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We calculate the two-loop effective action of QED for arbitrary constant elec- 



^ ■ Abstract 

o 

Q>^ ' tromagnetic fields at finite temperature T in the limit of T much smaller than the 

. electron mass. It is shown that in this regime the two-loop contribution always ex- 

ceeds the influence of the one-loop part due to the thermal excitation of the internal 
photon. As an application, we study light propagation and photon splitting in the 
presence of a magnetic background field at low temperature. We furthermore discover 
a thermally induced contribution to pair production in electric fields. 



1 Introduction 

In the low-energy sector of the theory, the effects of quantum electrodynamics can be 
summarized in an effective action, the Heisenberg-Euler action, which enlarges the classi- 
cal theory of electrodynamics by non-linear self-interactions of the electromagnetic field. 
Technically speaking, this effective action arises from integrating out the massive (high- 
energy) degrees of freedom of electrons and positrons. This program has successfully been 
carried out to two- loop order [|l], H, If. 

The inclusion of finite-temperature effects at the one-loop level has also been considered 
in various papers [|[ ^ H, 0, and the real-time [0] as well as the imaginary-time 

formalism finally arrived at congruent results. 

This paper is devoted to an investigation of the thermal QED effective action at the 
two-loop level. But contrary to the zero-temperature case, where the two-loop contribution 
only represents a 1%-correction to the one- loop effective action, we demonstrate that the 
thermal two-loop contribution is of a qualitatively different kind than the thermal one-loop 
part and exceeds the latter by far in the low-temperature domain. 

The simple physical reason for this is the following: at one loop, one takes only the 
massive electrons and positrons as virtual loop particles into account (cf. Fig. 0(a)). Due 



1 



to the mass gap in the fermion spectrum, a heat bath at temperatures much below the 
electron mass m can hardly excite higher fermion states. Hence, one expects thermal one- 
loop effects to be suppressed by the electron mass. In fact, in a low-temperature expansion 
of the thermal one-loop effective action [T^ , one finds that each term is accompanied by a 
factor of exp(— m/T), exhibiting an exponential damping for T — * 0. 





(b) 



Figure 1: Diagrammatic representation of the one- loop (a) and two- loop (b) contribution 
to the effective QED action. The fermionic double line represents the coupling to all orders 
to the external electromagnetic field. 

On the other hand, the two-loop contribution to the thermal effective action involves 
a virtual photon within the fermion loop (cf. Fig.|l|(b)). Since the photon is massless, a 
heat bath of arbitrarily low temperature can easily excite higher photon states, implying 
a comparably strong influence of thermal effects on the effective action. In Sec. 2, we are 
able to show that the dominant contribution to the thermal two-loop effective action in 
the low-temperature limit is proportional to T^/m^. This power-law behavior always wins 
out over the exponential damping of the one-loop case, leading to a two-loop dominance in 
the low-temperature domain. One might ask whether this inversion of the loop hierarchy 
signals the failure of perturbation theory for finite-temperature QED. But, of course, this 
is not the case, since the inclusion of a virtual photon does not "amplify" the two-loop 
graph and higher ones. Rather, calculating the one-loop graph should only be rated as an 
inconsistent truncation of the theory, since the one-loop approximation does not include all 
species of particles as virtual ones. Besides, effective field theory techniques indicate that 
the three-loop contribution is of the order of T^/m^ for T/m <^ 1, thereby obeying 
the usual loop hierarchy. 

The present paper is organized as follows: In Sec. 2, we present the calculation of 
the two-loop effective QED action at finite temperature employing the imaginary-time 
formalism and concentrating on the low-temperature limit. The outcome will be valid for 
slowly varying external fields of arbitrary strength. 

Section 3 is devoted to an investigation of light propagation at finite temperature. 
While, on the one hand, the well-known result for the velocity shift 6v ~ T^/m'^ is redis- 
covered ||12|, O, |l5i , we are also able to determine further contributions to the velocity 



shift arising from a non-trivial interplay between temperature and an additional magnetic 
background field. 

In Sec. 4, we study aspects of thermally induced photon splitting. Therein, we point 
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out that the thermal two-loop contribution to the splitting process exceeds the zero- 
temperature and one- loop contributions in the low-temperature and weak- field limit, but 
is negligible in comparison to other thermally induced scattering processes. 

Sections 3 and 4 are mainly concerned with the limit of a weak magnetic background 
field and low-frequency photons {uj <^ m), and therefore represent only a first glance at 
these extensive subjects. In fact, the quantitative results for this energy regime describe 
only tiny effects; a relevance for astrophysical topics such as pulsar physics has not been 
identified up to now. However, the intention of the present work is a more categoric 
one, namely, to elucidate the mechanism for a violation of the usual loop hierarchy of 
perturbative thermal field theories involving virtual massless particles. 

In Sec. 5, we calculate the thermal contribution to Schwinger's famous pair-production 
formula for constant electric background fields in the low-temperature limit. Here, a 
thermal one-loop contribution surprisingly does not exist 0, , since the thermal one- loop 
effective action is purely real by construction. Hence, the findings of Sec. 5 prove the 
existence of thermally induced pair production - an effect which has been searched for for 
15 years ^, |r^, |T^. In the low-temperature limit, we find that the situation of a 
strong electric field is dominated by the zero-temperature part (Schwinger formula), while 
the thermal contribution can become dominant for a weak electric field. Unfortunately, the 
experimentally more interesting high-temperature limit cannot be covered by our approach. 

One last word of caution: the inclusion of electric background fields in finite-temperature 
QED is always plagued with the question of how violently this collides with assumptions 
on thermal equilibrium. In fact, electric fields and thermal equilibrium exclude each other, 
thus questioning the physical meaning of the results of Sec. 5 at least quantitatively. How- 
ever, it is reasonable to assume the existence of an at least small window of parameters in 
the low-temperature and weak-field domain for which the thermal-equilibrium calculation 
represents a good approximation. Moreover, the knowledge of the effective Lagrangian 
including a full dependence on all possible field configurations is mandatory to derive 
equations of motion for the fields, even in the limit of vanishing electric fields. 



2 Two-Loop Effective Action of QED at Low Temper- 
ature 

In the following, we will outline the calculation of the two-loop effective action, concentrat- 
ing on the low-temperature limit where a two-loop dominance is expected. The calculation 
is necessarily very technical, wherefore some details are left for the appendices.^ 

But before we get down to business, it is useful to clarify our notation. From the field 
strength tensor F^*^ and its dual *Fi^p = \e^p^yF^^'' , we can construct the following standard 

^The primarily phenomenologically interested reader may just take notice of the following conventions 
(0)"®! then directly consult Eqs. (p9|)-(|3^), and skip the remainder of the section. 
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gauge and Lorentz invariants: 



g = ^F^'^*F^, = -E-B = ab, (1) 

where, for reasons of convenience, we also introduced the secular invariants 

a = \l + g2 5 = ^jr2 + g2 (2) 

and we assumed without loss of generality that a Lorentz system exists in which the electric 
and magnetic field are anti-parallel. In this particular frame, the secular invariants can be 
identified with the field strengths: a = B = |B|, b = E = |E|. 

When the physical system involves another vector, say, a momentum 4-vector fc^ = 
{k^, k), we can form another field invariant (metric: g = (— , +, +, +)): 

= |k|2s2sin2^B + Ikl^E^sin^^B- A;2e2 + 2A;°E- (kx B), (3) 

where 9b {Oe) denotes the angle between the magnetic (electric) field and the 3-space 
vector k. 

In relativistic equilibrium thermodynamics, temperature can be associated with the 
invariant norm of a 4-vector n^: n^n^ = — T^. On the other hand, is related to the 
4- velocity vector u'^ of the heat bath by: = T u^. E.g., in the heat-bath rest frame, 
takes the form: = (1, 0, 0, 0). Hence, we can introduce one further invariant (beside the 
temperature itself): 

£ = {u,Fn{u,F-'^). (4) 

E.g., in the heat-bath rest frame, S simply reduces to £^ = E^. Since the effective La- 
grangian is a Lorentz covariant and gauge-invariant quantity, it can only be a function 
of the complete set of invariants of the system under consideration. Hence, we expect a 
finite-temperature effective QED Lagrangian of the form: 

C = CiS,J^,g,T). (5) 

Equipped with these conventions, we now turn to the calculation. 

The two-loop contribution to the effective action/Lagrangian is generally given by 
the diagram in Fig. |l](b). This translates into the following formula in coordinate space 
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D^^{x-x'), (6) 

where G{x,x'\A) represents the fermionic Green's function for the Dirac operator in pres- 
ence of an external electromagnetic field A. D^i, denotes the photon propagator. Through- 
out the paper, we assume the background field to be constant or at least slowly varying 
compared to the scale of the Compton wavelength; therefore, the fermionic Green's function 
can be written as: 

G{x,x'\A) = I ^^e'p(^--'^g(p), (7) 

where g{p) denotes the Fourier transform of G{x, x'\A) depending only on the field strength, 
and ^{x, x') is the holonomy carrying the complete gauge dependence of the Green's func- 
tion. Inserting Eq. (|^ into Eq. (|^) leads us to the object x) = '^'(O)) where 
the right-hand side represents the holonomy evaluated for a closed path. For a simply 
connected manifold such as the Minkowski space, $(0) = 1; hence, it does not contribute 
to the zero-temperature Lagrangian. For a non-simply connected manifold such as the 
finite-temperature coordinate space (R x 5^), $(0) can pick up a winding number [§. 
However, in the present case, we restrict our considerations to a situation with zero density, 
which implies the existence of a gauge in which Aq = 0. Then, $(0) = 1 ^ind the influence 
of the holonomy can be discarded. 
This leads us to the representation 

for the two- loop Lagrangian, where D^^(k) denotes the photon propagator in momentum 
space, and we introduced the one-loop polarization tensor in an arbitrary constant external 
background field: 

W^{k) = -\e^ j ^^i,[Yg{p)^^g{p-k)\. (9) 

So we have finally arrived at the well-known fact that the two-loop effective action can 
be obtained from the polarization tensor in an external field by glueing the external lines 
together. 

The transition to finite-temperature field theory can now be made within the imaginary- 
time formalism by replacing the momentum integration over the zeroth component in Eqs. 
(^) and by a summation over bosonic and fermionic Matsubara frequencies, respectively. 



C^ = - i d'x'ti 



YG{x,x\A)YG{x,x\A) 
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E.g., performing this procedure in Eq. (^) corresponds to thermalizing the fermions in the 
loop. Now we come to an important point: confining ourselves to the low-temperature 
domain where T <C m, we know from the one-loop calculations [2C] that thermal 
fermionic effects are suppressed by factors of e"™'''^, indicating that the mass of the fermions 
suppresses thermal excitations. Hence, thermalizing the polarization tensor contributes at 
most terms of order e"*"/"^ to the two-loop Lagrangian for T <^ m; these are furthermore 
accompanied by an additional factor of the coupling constant a and can therefore be 
neglected compared to the one-loop terms. At low temperature, it is therefore sufficient to 
thermalize the internal photon only in order to obtain the leading T-dependence of C} . 
Since, in Feynman gauge, the photon propagator reads 



fc^ — ie ' 



+,+,+), 



(10) 



the introduction of bosonic Matsubara frequencies {k^Y 
us to:R 



C 



2+2T 



-(27rTn)2, n G Z, leads 



(27r)3 P - ie 



-Tl\{k). 



'Ill 



From now on, we write for the zero-temperature two-loop Lagrangian, C^^ for the purely 
thermal part, and £2+2T ^^j, their sum. In Eq. (|ll|), we need the trace of the polarization 
tensor in constant but otherwise arbitrary electromagnetic fields. In the literature, there 
are various equivalent representations for 11^,^. For the present purpose, it is useful to 
derive our own one which is based on a calculation of Urrutia ETI]. Details are presented 
in Appendix A. 

Inserting representation ( |A.19|) of the Appendix for II'^^ into Eq. ([Tl|) 
the Lagrangian: 



we obtain for 



C 



2+2T 



T a V ^ 

2" 2^^ 



oo 1 

(Pk f ds f du e~^^'^° eas ebs 



(27r)3 



2 a"^ + b'^ sin eas sinh ebs 



Zk 



fc2 



— (iV2 - iVi) + i2No{a' + b'') + b''N2 + a'iVi) 



(12) 



where the 0O) A^o? -Aj are functions of the integration variables s and v and of the invariants 
a and 6; only 0o depends additionally on Zk as defined in Eq. (^. Their explicit form can 

^Of course, the present calculation does not necessarily have to be performed in the imaginary-time 
formalism. E.g., instead of Eq. (p^), we could as well work with the real-time representation of the thermal 
photon propagator. We could even use the one-component formalism only, since we merely consider the 
photon to be thermalized. However, from our viewpoint, the calculations in the imaginary-time formalism 
appear a bit simpler since the momentum integrals will remain Gaussian. Of course, this might be just a 
matter of taste. 
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be looked up in Eqs. ( |A.16| ), ( |A.2CI| ) and ( [A.21| ). In order to ensure convergence of the 

ie for the mass term in 0o is 



m 



proper-time integrals, the causal prescription m 
understood; this agrees with deforming the s-contour slightly below the real axis. 

Now, the aim is to perform the /c-momentum integration/summation; note that the 
fc-dependence is contained in (po, Zk (and k"^, of course). Concentrating on this step, we 
encounter the integrals: 



h 



-iscpo 



(2^ 



(27r)3 F-ie 



(fc0)2 



(fc0)2 



(13) 



which allow us to write the Lagrangian (|T2D in terms of 



2+2T 



oo 1 

a f ds fdu 



eas ebs 



4vr 7 s J 2 (a^ + 6^) sineassinhefos 



{N2 - Ni) h + (2Aro(a' + &') + b''N2 + a^Ni)h 



(14) 



Employing Eq. ( [A.20|) for 0o, we can put down the evaluation of I2 to the one of Ii. 



d 



OA, 



(27r)3 P-ie 



dA, h, 



(15) 



where and Ak again are functions of the integration variables s and u and of the 
invariants a and b, and are defined in Eq. (|A.21| ). In view of Eq. (0), it is sufficient to 
consider the momentum integration/summation for Ji only: 



1 1 — Te 



E 



d^k 
(2^ 



(16) 



(fc0)2 



At this stage, the finite-temperature coordinate frame as introduced in ^ becomes ex- 
tremely useful, since it enables us to perform the calculation in terms of the invariants. 
This coordinate system is adapted to the situation of electromagnetic fields at finite tem- 
perature in a way that the components of any tensor-valued function of the field strength 
can be expressed in terms of the invariants S, JF, and Q. Again, details are presented in 
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the appendix (App. B), from where we take the final formula for the exponent of Eq. (|I6]) 
(cf. Eq. (m%): 



Ak + {a^-b^+S)A, 



(17) 



where k^, fc^, represent the components of the rotated momentum vector k"^ = e^f^k'^, 
and e^^ denotes the vierbein which mediates between the given coordinate system and the 
finite-temperature coordinate frame (cf. Eq. ( |B.23|) ). Since the transformation into the 
new reference frame is only a rigid rotation in Minkowski space, no Jacobian arises for the 
measure of the momentum integral. Hence, only integrals of Gaussian type are present in 
Eq. (|TB|), which can easily be performed to give: 



/i = T^^^^Ve-^-" (18) 

where it was convenient to introduce the short forms: 

qa:=Ak + a^A„ qb:=Ak-b^A„ p := Ak + {a^ -b^ +S)A,. (19) 
The sum in Eq. ([T^) can be rewritten with the aid of a Poisson resummation of the form: 

(20) 



exp(— a(ri — zY) = \ — exp ( — 2mzn j . 



With 2 = and a = we obtain for Eq. (H): 

-im?s 1 -\m?s i f 2 \ 

^^"^"" + ^' = W^ + V^^E-Pf-^^J' (21) 
IStt^ qaqh Svr^ g^g,, ^ V 4TV 

where we separated the (n = 0)-term from the remaining sum in order to find the (T = 0)- 
contribution. The first term in Eq. (^) [in = 0)-term), namely, is independent of T and 
S, while the second term vanishes in the limit T — > exponentially. In App. C, we check 
explicitly that the first term of Eq. ( pTD indeed leads to the (unrenormalized) two-loop 
Lagrangian for arbitrary constant electromagnetic fields at zero temperature. E.g., for 
purely magnetic fields, the representation of Dittrich and Renter M is rediscovered. 
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For our finite-temperature considerations, we will only keep the second term of Eq. 
(PI), which we denote by if in the following. Concerning the formula for C^^ in Eq. ([I^), 
Ii is already in its final form (it will turn out later that this term is subdominant in the 
low-T limit and only contains the important contributions). Hence, let us turn to the 



evaluation of , i.e., the thermal part of Eq. (|15]); for this, we have to interpret /f as a 



function of and Aj. (remember: Qa, Qb and p are functions of and A^): 

Z ^ Z J z 



Ak 



Upon a substitution of the integration variable,^ 



f22) 



where we defined the auxihary integrah 







u := '-^ (24) 

P [s' + qa){s' + qb) 

ds' du 



{s' + qa){s' + qb) + 2jMi^2p-qa-qb)u + {qa-qbfu^^ 

the auxiliary integral becomes: 

J{A-.) = f exp (-4,^ A ■ (25) 

''M + '-^{2p-qa-qb)u + {qa-qbyu^ V 4Tg,g, J 



Now we come to an important point: since we only thermalized the photons, our effective 
Lagrangian C"^^ is only valid for T <^ m anyway. Nevertheless, our formulas also contain 
information about the high-temperature domain which we should discard, since it is incom- 
plete. Regarding Eq. (pS]), the exponential function causes the integrand to be extremely 
small for small values of T, except where u is also small. Hence, the auxiliary integral is 
mainly determined by the lower end of the integration interval. 

Taking these considerations into account, we expand the square root for small values 
of u and then extend the integration interval to infinity (in fact, maintaining 1 as the 

^Resolving for s' — s'{u) leads to a quadratic equation from which the positive root has to be taken in 
order to take care of the integral boundaries. 
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upper bound only creates terms of the order exp(— (2nm)/T), which are subdominant in 
the low-temperature limit). The remaining w-integration can then easily be performed for 
each order in the w-expansion; up to m^, we obtain: 



(26) 



^ - J(^) = -25^(^ + S)- 2^^ ^g2_ 3(_^ _^ ^)2^)^^ ^ ^(TVnS). (27) 



Upon differentiation, the T^-dependence drops out, and we get (cf. Eq. (|T9| 



In this equation, we indeed discover a power-law dependence on the temperature, which 
will directly translate into a power-law dependence of the two-loop effective action after 
insertion into Eqs. (^) and ([T^) . Technically speaking, this arises from the fact that the 
omnipresent exponential factor exp{—^^-^^u), which finally causes exponential damping 
for T/m — > 0, becomes equal to 1 after the w-integration at the lower bound at u = 0. 

At this stage, it is important to observe that the ^-integration appears only in ij (via 
the ^'^-integration in Eq. ([13|) ) and not in if. Therefore, If will always contain exponential 
damping factors in the limit T ^ 0. Even the remaining proper-time integrations do not 
provide for a mechanism similar to the w-integration, since for large s, the mass factor 
exp(— im^s) with the causal prescription — > — ie causes the integrand to vanish, and 
for small s, the combination in the exponent becomes: 



QaQb 

P 4i 1 , , , 

qaQb S 

Obviously, inserting Eq. (pSj) into the exponent leads to an exponential fall off (bearing 
in mind that the s-contour will run slightly below the real axis). Similar conclusions can 
be drawn for the //-integration. To summarize these technical considerations, we conclude 
that only the term containing If (thermal part of I2) in Eq. (|14D contributes dominantly 
to C^^ in the low-temperature limit. 

Inserting the first and second term of j^J{Az) in Eq. (|27|) successively into Eq. ( ^2]) 



and then into Eq. (14), we obtain the dominant terms of order and of the two-loop 



effective QED Lagrangian at low temperature; particularly for the T^-term, different useful 
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representations can be given: 



2T 



eas ebs {N2 — N: 



00 1 

90 J s J 2 sin eas sinh ebs a? + 6^ 

-1 

00 



an 



-—T\J^ + S) / — 



45 



ds 1 



s a' 



+ 62 



1 — eas cot eas 1 — ehs coth e6s 

eos coth ehs h eas cot eas- 



sm eas 



a^ + b^ 

The term proportional to reads: 



sinh ebs 



00 

— e~'™'^'^eas cot eas ebs coth e6s 



in'^J s^ 




(29) 
.(30) 



2T 



945 



00 1 2 



s J 2 a^ + 52 gin eas sinh ebs 
-1 



(31) 



where Ni and are functions of the integration variables and the invariants a and b (not 
of S), and are defined in Eqs. ([A .161) and ( |A.21| ). The //-integration can be performed 
analytically, but the extensive result does not provide for new insights; hence we do not 
bother to write it down. 

These equations represent the central result of the present work; therefore, a few of 
their properties should be stressed: 

1) While we worked explicitly in the low-temperature approximation T <^ m, we put no 
restrictions on the strength of the electromagnetic fields. 

2) The low-temperature Lagrangians contain arbitrary powers of the invariants a and b 
(equivalently and Q), but the additional invariant at finite temperature S only appears 
linearly in the T'^-term and quadratic in the T^-term. The small-T expansion thus corre- 
sponds to a small-£^ expansion. 

3) The fact that only the integral /J with the prefactor (A'^2 ~ ^1) contributes to the low- 
temperature Lagrangian in Eq. {\L^) implies that only the spatially transversal modes Hn 
and of the polarization tensor ( |A.15| ) play a role in this thermalized virtual two-loop 
process. The time-like or longitudinal mode Ho (depending on the character of k^) might 
become important at higher values of temperature. 

4) The fact that the invariant S always appears in the combination J-' + S ensures a kind 
of dual invariance of the Lagrangian. Under the replacement E — B and B — E, the 
invariants change into JF — JF, Q ^ —Q and S ^ S + 2JF, so that J-' + S remains 
invariant. 
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5) The T*^-term of C^^ as exhibited in Eq. (|30| ) possesses the pecuharity of being derivable 
from the one-loop zero-temperature Lagrangian which we marked by square brackets in 
Eq. (|30D after the derivative terms. This will be elucidated a bit further in the following 
section. 

For the remainder of this section, we will discuss certain limiting cases of the two-loop 
low-temperature Lagrangian. First, let us concentrate on a weak-field expansion which 
corresponds to a small-s expansion of the proper-time integral due to the exponential 
mass factor. Expanding the integrands for small values of s (except the mass factor) and 
integrating over v and s, leads us to the dominant terms in the weak- field limit: 



^2T 



2^6 



2025 
213^3^5 ye 

36 ■ 5 ■ 72 ^ 



34 ■ 52 ■ 7 ^ 



0(3), 



■ (2J^' + ^ET + 3^2 _ (ry^^^^ 

mr 



(32) 
(33) 



where 0(3) signals that we omitted terms of third order in the field invariants (sixth order 
in the field strength). Note that no linear term in the field invariants to order exists. For 
the terms of quadratic order, the T^-term is subdominant for T /m < 0.05, and amounts up 
to a 10%-correction to the T^-term for T/m ~ 0.1. For even larger values of temperature, 
we expect the failure of the low-temperature approximation. 

Finally, we consider C'^'^\rp4, in the limit of purely magnetic background fields: 6^0, 
a ^ B , + £ ^ ^tB"^ . The T^-term in Eq. (pOj) then reduces to: 



c:''\B) 



rp4 



an 
90 



rp4 



— e ' 

z 



1 — z coth z 1 

2 1 — zcothz 

sinh z 3 



(34) 



where we have performed the substitution eas — > —iz in concordance with the causal 



prescription m 



m 



ie. Incidentally, the limit of purely electric fields can simply be 



obtained by replacing B —lE and multiplying Eq. (0) by (—1). 

Introducing the critical field strength := we can evaluate the integral in Eq. 
(P^) analytically |jl5[Q and obtain: 



c:''{B) 



^4 



90 



2S^ 3 



i)V^(l + |-)-^lnr(|-) 



35| 

452 



.:^ + :^ln27r+-+4C'(-l 
2B B 6 ^ ^ 



BcT ) J 

' ^ SB, 



, (35) 



where iIj{x) denotes the logarithmic derivative of the F-function, and C'{s,q) is the first 
derivative of the Hurwitz (^-function with respect to its first argument. 

^We take the opportunity to remark that there is a misprint in the corresponding integration result in 
Ipf; the term (+1/3) has to be replaced by (+1/6) (cf. Eq. (|3l 
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For strong magnetic fields, B ^ Bcr, the last term in square brackets in Eq. ( pS] ) 
dominates the whole expression, and we find a linear increase of the effective Lagrangian: 



rp4 



an ^4 eB 



270 



(36) 



This contribution remains subdominant compared to the one arising from pure vacuum 
polarization ~ B^ln^, which is not astonishing, since the magnetization of (real) ther- 
malized plasma particles is bounded: the spins can maximally be completely aligned. In 
contrast, the non-linearities of vacuum polarization set no such upper bound. Quantita- 



tively, the same result was found for the thermal one-loop contribution [53 . 



3 Light Propagation 

As a first application, we study the propagation of plane light waves at finite temperature 
and in a magnetic background. The subject of light propagation has recently gained 
renewed interest due to its accessibility to current experimental facilities . 

In the limit of light of low-frequency -C m, the effective action for slowly varying fields 
has proved useful for obtaining velocity shifts, i.e., refractive indices of QED vacua which 



are modified by various external perturbations such as fields and temperature [0, ^ 
In this limit of low frequencies and smooth external perturbations, the terms involving 
derivatives of the fields in a derivative expansion of the effective action can be neglected, 
and the constant-field approximation is appropriate. 



The case of light propagation at finite temperature has been investigated in |20[ from 
a general viewpoint for a class of Lagrangians depending on the invariants S, J-", Q,T in an 
arbitrary way. Therein, a light cone condition representing a sum-rule for the polarization 
modes of the propagating light has been derived; this has been exploited for a detailed 
investigation of light propagation at finite temperature to one-loop order by an insertion of 
the thermal one-loop effective Lagrangian of QED. It has been emphasized that these one- 
loop studies apply to a domain of intermediate values of temperature ~ 0.1 < T/m <~ 1, 
where two-loop as well as plasma effects remain subdominant. 



The famous results for the low-temperature velocity shift 5v ~ T/m |T2|, 
could not have been rediscovered by this first-principle investigation, because the thermal 
two-loop effective action was not at hand. In the present work, we intend to fill this last 
gap. 

Let us first consider the situation of a thermalized QED vacuum without an additional 
background field. In the low-temperature domain, this vacuum is then characterized by 
the Lagrangian C = —J-'+C'^^, where — represents the classical Maxwell term. Following 
the lines of the phase and group velocity of a propagating plane wave is then given 
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by: 

= . I \d,c > (37) 

where f = is constructed from the wave vector of the propagating hght, and it is 
understood that the partial derivatives of C are evaluated in the zero-field limit. Inserting 
Eqs. (H) and (|3|) into Eq. (13), leads us to: 

^ 2025 m4 

Note that there is no T^-term, since C?'^\tg is at least quadratic in the field invariants. In 
Eq. (|38[) , we rediscovered the well-known velocity shifts for light propagation in a thermal 
background as found in []T^ |T3| via the two-loop polarization operator and in via 



considering vacuum expectation values of field bilinears in a thermal background. The 
here-presented rederivation within the effective action approach from first principles thus 
can be viewed as an independent check of our calculations and of the light cone condition 
as derived in |j20|. 



But we can go one step further and additionally take a weak external magnetic field 
into account; the light cone condition in this case reads 



0= {d:FC-d£C-J^dlC)k^+]^{dl+dl)Czk + 2d£C{kuf, (39) 

where denotes the 4- velocity vector of the heat bath and Zk is defined in Eq. (^. The 
Lagrangian describing a thermal QED vacuum with weak magnetic background fields at 
finite temperature is given hy C = —T^C} + where C} denotes the one-loop effective 
Lagrangian at zero temperature. Up to the second order in the invariants, this famous 
Heisenberg-Euler Lagrangian C} is given by: 

= + -g^. 40 

45 m4 45 m^^ ^ ^ 



Inserting all the relevant contributions to C into the light cone condition Eq. (pOf) , the 
light velocity to lowest order in the parameters T and B finally yields: 

v^= 1 jB^sin^ 9b -2 aV— + {— "tt^^ 5^ 1 + sin^ , 

45 m4 ^ 2025 45 ^3^ ■ 5 ■ 7 ■ 11 my ^ ^' 

(41) 

where 6b denotes the angle between the propagation direction and the magnetic field (cf. 
Eq. (PP). The second and third term are the well-known velocity shifts for purely magnetic 
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25| , p6| and purely thermal vacua (cf. Eq. (^)), respectively. The last term describes 



a non-trivial interplay between these two vacuum modifications. The latter can best be 
elucidated in the various limits of the angle Ob'-, for orthogonal propagation to the magnetic 
field 6b = vr/2, we get: 

= 1- 2 aV— -B^ 1 - 0.15... ■—:]. 42 

2025 45m4 V ^ ' ) ^ ' 

For parallel propagation to the magnetic field = 0, we find: 



(43) 




2025 m4 



Since T/m and eB /m? are considered to be small in each case, the corrections to the pure 
effects in the mixed situation are comparably small. Note that the mixed thermal and 
magnetic corrections always diminish the values for the velocity shift of the pure magnetic 
or thermal situations. Let us finally remind the reader that the here-given velocities hold for 
low-frequency light {uj ^ m) only, and represent averages over the two possible polarization 
modes. While for the purely thermal case the polarization modes cannot be distinguished, 
the situation involving an electromagnetic field generally leads to birefringence due to the 
existence of a preferred direction of the field lines. 



Let us finally comment on the earlier works [jT3|, |T5[ related to the issue of light prop- 
agation in a thermal background. The philosophy therein was to calculate the velocity 
shifts in a purely (weak) electromagnetic background first, and then take thermal vacuum 
expectation values of the field bilinears. Expressing this in formulas, we first recall the 
expression for the propagation-direction-averaged light velocity in a weak electromagnetic 



background from |15|: 



V 



3 



where = \{E'^ + B"^) denotes the 00-component of the energy-momentum tensor, i.e., 
energy density of the electromagnetic field. In the weak- field limit, ((9|r + d'g)C is field 
independent: 2||^ (cf. Eq. (0)); therefore, taking thermal vacuum expectation values 



of the field quantities in Eq. (^J) is simply equivalent to replacing T by (T y = j^T ■ 
This then leads to the correct result as given in Eq. (B3). 



From the viewpoint of the present work, the correctness of the approach of [|T3|, [15 1 
arises from the special form of the low-temperature two-loop Lagrangian C'^'^It^, as given 
in Eq. (|30D. Since 

^2(91 + d!) = d'^ + dl (45) 
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Eq. (pOD can also be written as: 



2T 



rp4 



(46) 



Incidentally, Eq. (^) holds for arbitrary field strength, but, in this line of argument, it is 
required for weak fields only. Inserting Eq. (^6]) into the correct light cone condition at 
finite temperature, i.e., Eq. (|37D, we obtain to lowest order: 



1 - 2drC = 1 



/rrOO\T 



y 



(47) 



which is equal to the heuristically deduced light cone condition for a thermal QED vacuum 

in, Hi. 

Note that the combined low-temperature/weak-field effects as given in Eqs. (^T])-(^3|) 
could not have been found in [|l5l, since the invariant structure is not completely taken 



into account in the heuristic approach. Whether the investigation of the intermediate- 
temperature domain to two-loop has been correctly modeled with the heuristic approach 
cannot be judged within the present work. Note, however, that the intermediate- 



m 



temperature domain is controlled by one-loop effects, leading to a maximum velocity shift 
of — 5fmax = As has been shown therein, the two-loop dominance is lost for 

T/m > 0.058. 



4 Photon Splitting 

Photon splitting in magnetic fields at zero temperature has been discussed comprehensively 
by Adler |2^, stressing its relevance for the photon physics of compact astrophysical objects 
(see also ||2^). For the description of the splitting process for low- frequency photons with 
uj <^ m at weak magnetic fields ^ <C 1, the use of the one- loop effective Lagrangian 
for weak fields is sufficient for obtaining a good estimate of the absorption coefficient for 
photon splitting. To be precise, the lowest order contribution to the splitting process comes 
from the terms of third order in the invariants (sixth order in the field strength) of C^, i.e., 
the hexagon graph with one incoming, two outgoing photons and three couplings to the 
external magnetic field. Neglecting dispersion effects, the box graph vanishes because of 
depending on JF and Q only, and because of the Lorentz kinematics of the photons.0 
The question of thermally induced photon splitting has recently been investigated by 
Elmfors and Skagerstam [|1^ with the aid of the thermal one-loop effective QED La- 
grangian; their studies were motivated by the fact that a vacuum may be a bad approx- 
imation for the surroundings of some astrophysical compact objects, while a thermalized 
environment at zero or finite density might be more appropriate. It turned out that, at 

''Taking dispersion effects into account, the box graph stiU is only an order a correction to the hexagon 
graph. 
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temperatures and magnetic fields at the scale of the electron mass, the thermal contri- 
bution can exceed the zero-temperature one, but these effects then are superimposed by 
Compton scattering of the photons with the plasma. In realistic situations, the thermally 
induced process will thus be of subdominant importance. 

In the following, we intend to complete these results about thermally induced photon 
splitting with the dominant low-temperature contributions stemming from the two-loop 
process. Hereby, we also concentrate on the splitting process (± ||i + 112)5 where a 
photon, with its electric field vector orthogonal (_L) to the plane spanned by the external 
magnetic field and the propagation direction, splits into two photons with their electric 
field vectors within (||) that plane.0 This is the only allowed process when dispersion effects 
are taken into account. 

As pointed out in [|10], the box-graph no longer vanishes at finite temperature, since 
the Lagrangian now involves an additional invariant. Hence, the lowest-order contribution 
to the photon-splitting matrix element is already produced by the terms of quadratic order 
in the invariants in Eqs. (|32D and ([33|). 

Without going into details, we recall that the splitting amplitude is obtained by attach- 
ing the external photon legs to the fermion loop, i.e., differentiating the effective action 
(which is represented by the loop) thrice with respect to the fields and then contracting 
the result with the field strengths of the involved photons. Hereby, one has to take into 
account that the effective Lagrangian now depends on three field invariants: S,J^, and Q. 
The thermal amplitude arising from the box-graph finally yields: 

A^(± ^ 111 + II2) = 2ujuJiUJ2Bsm9Bd£^C, (48) 

where lu,ui,uj2 denote the frequencies of the incoming and the two outgoing photons, 
respectively, and 9b again represents the angle between the propagation direction and the 
magnetic field. From the splitting amplitude, we obtain the absorption coefficient k via 
the formula: 



K = 

3271UJ 



— - / duji / duj2S(uj - ui - UJ2) M'^. (49) 
J J 



Inserting Eq. (|45|) for the thermal splitting amplitude into Eq. (H3) leads us to 



m 



2^ • 3 • Svr^ \ m? J \m/ 



sin^^B - {de:pCYm\ (50) 



Here, we encounter the typical (cu/m) ^-dependence of the photon-splitting absorption co- 
efficient for low-frequency photons. The appearance of the magnetic field to the second 



^Note that Adler's definition for the ||,_L-mode rely on the direction of the magnetic field vector of the 
photon and thus are opposite to ours. 



17 



power is directly related to the fact that the box-graph exhibits only one coupling to the 
external field. In contrast, Adler's result for the absorption coefficient at zero temperature 
arising from the hexagon graph reads |25]: 



T=0 



a 



m 



35 . 53 . 72 ^2 



eB 



-) 



(51) 



Here, the three couplings to the external magnetic field produce a B^- dependence of the 
absorption coefficient. Therefore, any finite-temperature contribution will exceed the zero- 
temperature one for small enough magnetic fields; but, of course, the absorption coefficients 
may then become very tiny. 

In order to obtain the one-loop and two-loop absorption coefficients for thermally in- 
duced photon splitting at low temperature, the derivatives of the corresponding Lagrangian 
are required in Eq. (|50D: 



"45" 



+ 



34-52 
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(52) 
(53) 



where we made use of the results of [0 for the low-temperature/weak- field approximation 
of the one-loop Lagrangian C^^ , and employed Eqs. ( p2D and (^) for the two-loop one. 
Obviously, inserting the two-loop terms from Eq. ( ^31) into Eq. (^) leads to a power-law 
dependence of the absorption coefficient ~ T^/m^, while the one-loop terms from Eq. (^2]) 
imply an exponential mass damping exp(— 2m/T) for T —>■ 0. 

As mentioned above, photons of frequencies below the pair-production threshold are 
not only exposed to splitting at finite temperature, but can also scatter directly with 
the plasma of electrons and positrons. Following [|T^], the absorption coefficient for the 
Compton process is given by 



m m TT^ 





(54) 



where oOg denotes the fermion energy Ue = ^/ p'^ + m^, and the cross section cxc for unpo- 
larized photons at tu/m ~ 1 is approximately given by: 



Ana' 
3m2 



(55) 



Although uj/m 1 formally represents the maximal limit of validity of our constant- 
field approximation for the effective action, we will continue to consider photons of that 
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frequency in the following, since, on the one hand, this circumvents a suppression of the 
absorption coefficients by the common factor {u/m)^, and on the other hand, it has been 



shown for the hexagon graph in |25| that the difference between u/m = 1- and uj/m ~ 0- 
calculations is negligible for weak magnetic fields. 

Finally, we have to consider another scattering process which arises from the presence 
of a heat bath: photon-photon scattering between the propagating photon and the black- 
body radiation of the thermal background. We estimate the absorption coefficient for this 
process by 

m m ' 

where denotes the density of photons and is given by: 

(57) 

Here we encounter the Riemannian ^-function with C(3) ~ 1.202. The total polarization- 
averaged cross-section for photon-photon scattering at low frequencies, as one obtains, e.g., 
from the Heisenberg-Euler Lagrangian |]28|, reads: 



_ 973 (^cm\^ 
~ 10125 \ 



(58) 



where cucm denotes the frequency of both photons in the center-of-mass frame. In order 
to determine cucm, we first have to find the mean frequency at temperature T. Averaging 
over the thermal probability distribution for the photons, we find the mean value lot = 

4 

30C(3) '^ ~ 2.701T. According to relativistic kinematics, the average value for the CM- 



frequency ujqm is given by ujqm = ^/ujujt/2 ~ 1.16\/Ta;, where we averaged over the 
propagation direction of the thermal photons. Putting everything together, we obtain for 
the absorption coefficient for photon-photon scattering with the thermal background: 

--rr ^ J-f\^^' (L)' ^ 5.21 • IQ-" f (^Y . (59) 

m 2^ ■ 3' ■ 5° ('(3)"' \m J \m/ \m J \m/ 

Since the average frequency of the heat-bath photons is proportional to the temperature, 
this formula becomes invalid for T ~ m and above, because we employed the low-frequency 
cross-section in Eq. (0). 

It is already clear from a qualitative viewpoint that there must be a domain where the 
two-loop splitting process at least exceeds the one-loop and the Compton contributions 
due to the power-law dependence on the temperature. But since k^'^ ~ (T/mY and 
~ (T/mY, the two-loop contribution will eventually be surpassed by the photon- 
photon scattering for T 0. 
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(b): ^ = 10-^ ^ 
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Figure 2: Absorption coefficient /t in units of the electron mass versus temperature T in 
units of the electron mass. In Fig. (a), the various contributions are plotted for parameter 
values of a realistic astrophysical system. In Fig. (b), the parameters are chosen in such a 
way that the two-loop dominance over the one-loop and the Compton process is revealed; 
the photon-photon scattering contribution cannot be overtaken in the low-temperure limit. 

However, quantitative results can only be revealed by numerical studies. In fact, as 
shown in Fig. |^(a), the two- loop contribution is completely irrelevant for parameter values 
which may be appropriate for a neutron star system and which are close to the upper bound 
of vahdity of our approximation: ^ = 0.2, u/m = 1, sin6'B = 1, and T/m = 0.05 . . . 0.1. 
Even the one-loop contribution is small compared to the zero-temperature result; but all 
are negligible compared to the Compton process. 

Concentrating on the relative strengths of the thermal splitting processes, the one- 
loop contribution loses its major role for T/m < 0.041, where its exponential decrease is 
surpassed by the two-loop power law. 

In order to find a domain in which the two-loop splitting wins out over the zero- 
temperature process, we have to look at smaller values of the magnetic field strength; e.g., 
at values of temperature T/m = 0.025, the two- loop process exceeds the zero-temperature 
one for ^ < 2.1 ■ 10~^. Since these are more moderate field strengths, the absorption 
coefficient naturally becomes very small: n/m ~ 10~^^ . . . 10^^^. Hence, in order to be 
able to measure the splitting rate, the extension of the magnetic field in which the photon 
propagates must be comparable to galactic scales. 

Finally, we have plotted the Compton and photon-photon absorption coefficients, hq 
and K^:y, and the two-loop coefficient k"^^ for a weak magnetic field ^ = 10~^ at T/m = 
0.001... 0.1 in Fig. 0(b). Obviously, the Compton process loses its dominant role for 
T/m < 0.03; below, the absorption coefficient is ruled by the photon-photon scattering as 
long as the temperature does not become so small that only the zero-temperature amplitude 
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remains. As is also made visible in Fig. 0(b), the two- loop contribution does not exceed the 
photon-photon process, due to the weaker temperature dependence of the latter. Hence, 
we may summarize that the photon absorption coefficient in the low-temperature domain 
is either dominated by the zero-temperature contribution for strong magnetic fields or by 
the photon-photon scattering with the thermal background for weak fields. So the two-loop 
contribution always belongs to the top fiight but is never ranked first. 

In order to account for realistic astrophysical systems, it is compulsory to include a 
finite chemical potential. First estimates can be found in to one-loop order, where 
signals have been found that a finite chemical potential of /x ~ m may induce an increase 
of the thermal splitting amplitude at low temperatures. In order to settle this question 
properly, the present paper shows that an investigation of these systems should take the 
two- loop contributions into account. 

Let us conclude this section with the remark that in order to obtain the sum of the zero- 
temperature and the thermal contributions to the photon splitting absorption coefficient, 
the amplitudes must be summed up coherently, since the final states of the processes 
coincide, and the thermal vacuum with a constant background field does not provide for a 
mechanism of decoherence. While the zero-temperature amplitude as well as the thermal 
one- loop amplitude are strictly positive, the T'^-term in Eq. (53) contributes with a negative 
sign. Hence, an exceptional curve in the parameter space of ^ and T/m exists where the 
thermal two-loop amplitude interferes with the thermal one-loop and zero-temperature 
amplitudes destructively so that photon splitting vanishes. 



5 Pair Production 

Thermally induced pair production in electric fields has been searched for at the one-loop 
level for a long time H, |, |l^ |18|, |l9l with extremely contrary results. In our opinion, 
the final concordant judgement in the real-time formalism [0], the functional Schrodinger 
approach [T^, as well as the imaginary-time formalism is that there is no imaginary 



part in the thermal contribution to the effective action to one loop, implying the absence of 
thermally induced pair production to this order of calculation. As already mentioned in the 
introduction, drawing conclusions from an imaginary part of the thermal effective action 
to pair production is not as immediate and straightforward as at zero-temperature, since 
the presence of an electric pair-producing field and the thermal equilibrium assumption 
which is inherent to our approach contradict each other. 

In the following, we simply assume that on the one hand, the time scale of pair pro- 
duction is much shorter than the time scale of the depletion of the electric field so that 
dynamical back-reactions can be neglected (this assumption is familiar from the zero- 
temperature Schwinger formula). On the other hand, we also assume that the state of the 
plasma is appropriately approximated by a thermal equilibrium although it is exposed to 
an electric field. Whether the assumption on thermal equilibrium is justified in concrete 
experimental situations such as, e.g., heavy ion collisions, is still under discussion. 

Let us now turn to the computation of the imaginary part of the two-loop thermal 
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effective action for external electric fields. For this, we concentrate on the T'^-contribution 
as given in Eq. (^91). For purely electric fields, a^O,b-^E,S + J^-^ this reads: 



^4 



90 



1 1 — z coth z 

-z coth z H 2 

3 sinh z 



(60) 



where we substituted z = eEs. For reasons of convenience, it is useful to abbreviate 



V 



:= which denotes the dimensionless ratio between the electric field and the critical 



field strength E, 
C^iE) 



Integrating the 1 / sinh^ ^-term by parts leads us to: 



T lim< 



2^4 



90 



-0 1 2e2 2 4r/2 



2 + + —2+ I ^''~'H-S---^ + ^2)^^'^'y (61) 



Here, it should be pointed out that the isolated pole in the first term of the curly brackets 
does not signal a divergence, but simply cancels the pole at the lower bound of the integral; 
the whole expression is still finite. Our aim is to evaluate the imaginary part of Eq. ([61|); 
for this, the behavior of the integral at the lower bound is of no interest. An imaginary part 
Im C'^'^ {E)\t4, arises from the poles of the coth2;-term on the imaginary axis aX z = iiravr, 
n = 1,2,.... 

Decomposing the exponential function into cos +i sin, it becomes obvious that the imag- 
inary parts of the integrand are even functions in z, while the real parts are odd. Thus, 
extending the integration interval from — oo to oo exactly cancels the real parts and simply 
doubles the imaginary parts. We finally get: 
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(62) 



Now we can close the contour in the lower complex half plane, which is in agreement with 



the causal prescription m 



ie. The value of the integral is then simply given by 



the sum of the residues of the cothz-poles at 2; = — ivm, n = 1, 2 
ImC^iE) 



Hence, we arrive at: 



ttTT^ . ^ / 1 1 1 1 \ eE 

^-^ \6 nrcr] n"^7r^ 2r]^y 
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n=l 



which represents our final result for the imaginary part of the thermal effective QED action 
at low temperature, and should be read side by side with Schwinger's one-loop result: 



Im C\E) 
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(64) 
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The sum in Eqs. (|^) and ( [6^ ) can be carried out analytically; but here, it should be 
sufficient to consider the limiting cases of weak and strong electric fields. 

In the weak-field limit, i.e., for small values of rj, the sum over n in Eq. (|63|) is dominated 
by the first term n = 1. Furthermore, it is the last term which is the most important one 
in parentheses. These considerations then lead us to: 

ImC'^ieE « m^) ^ ^ (65) 

180 1]^ 

Combining this with the weak-field approximation of Eq. (0), we get roughly for the total 
imaginary part of the effective Lagrangian: 

lmC{eE « m^) = mV^/^ (^ + ^ mV^/^ (4 ■ lO'V + 4 ■ IQ-^'^^^] . 

\87r'^ 180 1]'' m^J \ V J 

(66) 

E.g., for TjvTL':^ 0.1, where the present low-temperature approximation should still be 
appropriate, the thermal contribution can be neglected for ri > 0.1; both contributions 
become roughly equal for rj ^ 0.056 (and T/m = 0.1). For weaker fields and T/m ~ 0.1, 
the thermal contribution even becomes the dominant one. 

In the opposite limit, where r] ^ 1, i.e., for strong electric fields beyond the critical 
field strength, the 1/3 in parentheses dominates the expression in Eq. (|63|) , which then 
gives: 

ImC'^ieE » m') = ^T^V^ fe'^M" = = ^T'r^ + 0{rf). (67) 

^ ^ 270 ^ V / 270 1 - 270 ' w y v ; 

71=1 

Together with the strong-field approximation of the Schwinger formula, this gives: 

Im C{eE > m^) = r/ f — + — ^ r/ (^6.6 ■ IQ^^r] + 8.5 ■ 10"^^^ . (68) 



^487r 270 mV V 

Since Eq. (^) is valid for 77 ^ 1 and T/m <^ 1, the low-temperature contribution to 
\m.C{E) can be neglected for strong electric fields. Similarly to the case of strong magnetic 
fields, we find that the non-linearities of pure (zero-T) vacuum polarization exceed the 
polarizability of the thermally induced real plasma by far in the strong field limit. 

Nevertheless, in the limit of weak electric fields, thermal effects can increase the pair- 
production probability P = 1 — exp(— 2Im£(ii^)) significantly, as was shown in Eq. (p6|). 
Of course, for these values of 77, the total imaginary part is very small due to the inverse 
power of 7] in the exponential. 

Since we did not consider thermalized fermions, our approach is not capable of describ- 
ing high-temperature pair production, which would be desirable for forthcoming heavy-ion 
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collision experiments. However, as can be read off from our results for light propagation 
and photon splitting, extrapolating the power-law behavior to higher temperature scales of 
T ~ m or even T/m 3> 1 overestimates a possible two-loop contribution by far, since, for 
these values of temperature, the one-loop contribution can be expected to be the dominant 
one. The latter increases at most logarithmically with T. 

Therefore, it is reasonable to assume that the pair-production probability also increases 
at most logarithmically with T. In view of these considerations, a power-law growth as 



suggested in R 18, 13] does not appear plausible. Of course, in order to decide this question. 



the two-loop calculation has to be carried out for arbitrary values of temperature. 



6 Discussion 

In the present work, we calculated the thermal two-loop contribution to the effective QED 
action for arbitrary constant electromagnetic fields in the low-temperature limit, T/m ^ 1. 
Contrary to the usual loop hierarchy in a perturbation theory with small coupling, the 
thermal two-loop part is found to be dominating over the thermal one-loop part in the 
low-temperature limit, since the former exhibits a power-law behavior in T/m, while the 
latter is exponentially suppressed by factors of exp(— m/T). The physical reason behind 
this is that the one-loop approximation does not involve virtual photons, which, due to 
their being massless, can be more easily excited at low temperatures than massive fermions; 
thus, the one-loop approximation should be rated as an inconsistent truncation of finite- 
temperature QED for T much below the electron mass. 

The power-law dependence of the thermal effective action to two loop starting with 
T^/m^ implies a two-loop dominance in the low-energy domain of thermal QED, which 
holds roughly up to T/m ~ 0.05. 

For the subject of light propagation at finite temperature, this two-loop dominance has 
been known for some time from studies of the polarization tensor [p!2| , |T^. Moreover, for 
the subject of QED in a Casimir vacuum like the parallel-plate configuration, the two- 
loop dominance is very natural and well known, since the fermions are not considered to 
be subject to the periodic boundary conditions anyway. This gives rise to a non-trivial 
check of our results, since Casimir and finite-temperature calculations highly resemble each 
other. Replacing, as usual, T by l/(2a) in Eq. (|65D for the weak-field limit of the imaginary 
part of the effective Lagrangian, where a denotes the separation of the Casimir plates, we 
obtain: 



Im£'"(E) 



\ eE > ' ^ ' 



which agrees precisely with the findings of for the Casimir corrections to the Schwinger 
formula^]. 

^Actually, Eq. (|6i) agrees with the findings of except for a global sign; however, as was pointed 
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In order to illustrate the two-loop dominance, we studied light propagation and photon 
splitting in a weak magnetic background at low temperature. Since we are dealing with 
the two-loop level, the here-considered effects are naturally very tiny and a significant 
influence on, e.g., photon physics near astrophysical compact objects appears not very 
probable. One should rather take a closer look at photon physics on large galatic scales. 

Furthermore, we calculated the imaginary part of the thermal two-loop effective action 
for electric background fields at low temperature. Under mild assumptions, this result 
can be related to a thermally induced production probability of electron-positron pairs. 
Especially in the weak-field limit, the thermal contribution has a significant influence on 
the production rate. Since no thermal one-loop imaginary part exists, any finite two-loop 
result automatically dominates at any temperature scale. 

For the subjects of light propagation and photon splitting, the loop hierarchy is restored 
above T/m ~ 0.05. Already at this comparably low value of temperature, the thermal 
excitation of the fermions begins to compete with that of the virtual photon. Hence, a 
calculation of the two-loop thermal Lagrangian at intermediate or high temperatures would 
appear as an imposition, were it not for the high-temperature pair-production probability 
which is beyond the range of the one- loop approximation and of great interest for, e.g., 
heavy-ion collisions. 



Appendix 

A One-loop Polarization Tensor 

While the polarization tensor in an external magnetic field has been considered by many 
authors (a comprehensive study can, e.g., be found in |^), a generahzation to arbitrary 
constant electromagnetic fields in a straightforward manner is associated with a substantial 
increase in calculational difficulties. The problem was first solved by Batalin and Shabad 
3^ ; their extensive result was later brought into a practical form by Artimovich ^3 



In the following, we will briefiy sketch a simpler derivation of the polarization tensor in 
arbitrary constant electromagnetic fields; our approach is based on the findings of Urrutia 
Il2l| , who solved the problem for the special case of parallel electric and magnetic fields. 



Assume that the (— i?)-field and the B field point along the 3-axis. 4-vectors like the 
external momentum (cf. Fig. can then be decomposed into 

F = A;f + A;^, A;f = (fc°,0,0,fc^), k"^ = {0, k\ k\0). (A.l) 

In the same manner, tensors can be decomposed, e.g., g^'^ = gi^'^ + g^ . With respect to 

out by one of the authors in a footnote of the expression in |^ is wrong by a minus sign, which saves 
the day. 
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Figure 3: Diagrammatic representation of the one-loop polarization tensor. The fermionic 
double line represents the coupling to all orders to the external electromagnetic field. 



each subspace, we easily find the unique orthogonal vector to a given one: 



(A.2) 



Following Urrutia [21|, the polarization tensor for the special field configuration can be 
written as: 



oo 1 

^ ' ^ 2tt s 2 ] sinzsinh^' 



-1 



+ {g^^kl - k^lkl)N2 - {k^J^l + k^kDN^ 



+ c.t. \. (A.3) 



The electric and magnetic field strengths E, B are contained in the variables z := eBs and 
z' := eEs. The exponent 0o is given 



m + 



^if cosh z' — cosh vz' fc^ cos vz — cos z 



z' sinh z' 



+ 



z sm z 



The functions iVj read:[| 



(A.4) 



Nq = cosh uz' cos uz — sinh uz' sin uz cot z coth z' , 

=: iVi - iVo, 

=: N2 - No, 



cosh z' — cosh z/z' 
iVi = 2cosz : — 5 No 



sinh 



^ , , cos z/z — cos z 
No = 2 cosh 2' ^ N( 



sin^ z 







1 — cos z cos uz cosh z/2;' cosh z' — 1 , 
Ns = : — — : h sin uz sinh uz , 



sm z 



sinh 2;' 



(A.5) 



^This formula has been misprinted in Ref. |2T[ . 

^iVa differs from Urrutia's findings by a minus sign, since he considers parallel E- and i3-fields. 
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where we have incidentally defined the functions iVi 2 for later use. The determination of 
the contact term corresponds to a charge and field strength renormalization and yields: 

c.t. = -e-''"''(l - [g^'^e - k^'k") . (A.6) 



Now, one can show |Q that the Lorentz invariant form of the polarization tensor for 
arbitrary constant electromagnetic fields can be completely reconstructed from the special 
form given above for anti-parallel electric and magnetic fields. This is achieved by, first, a 
one-to-one mapping between Urrutia's scalar variables (fcy , /c^, B) and a set of invariants 
which reduce to Urrutia's variables in the special system: 

a ^ B, b ^ -E, 

Zk ^ -E'kl + B^kl, k'^kl + kl (A.7) 

The inverse map is obtained by a simple calculation; the non-trivial relations are: 

Secondly, the reconstruction requires a one-to-one mapping between Urrutia's tensor struc- 
tures in Eq. (|A.3|) and Lorentz covariant tensors which reduce to Urrutia's in the special 



system. For this, we need to introduce the following definitions. First, we employ a set of 
four linearly independent 4-vectors: 

k^", Fk^ = F^'''k^, F^k^" = F^'^'F^pk'^, and *Fk^ = *F^"ka. (A.9) 

From these, we construct the 4-vectors: 

< := , ^ ,J a*Fk^'-bFk^') k^;, 
" + b'' \ / " 

< := —^(b^Fk^ + aFk^") k'f , (A.IO) 

+ b'' \ / 

where the subscripts || and _L are to remind us of the meaning of v\\ and v± in the spe- 
cial Lorentz system (longitudinal and transversal part of k). Incidentally, they obey the 
relations (cf. Eq. ( |A.8| )): 



2 u a^k"^ Zk 2 u b'^k'^ -\- Zk „ /a ii\ 

^11 = Vll^ = ~ q2 ^ 52 ' ^± = = ^2 ^ 52 ' ^f^±M = 0. (A.ll) 
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Finally introducing the projectors 

1 



A;2 



P; 



„,2 ' 



+ g^- (n) 



,,2 ' 



which satisfy the usual projector identities, P^^ 
the one-to-one mapping: 



0,11, 



P 



0,11 



^o,||,±% 



(A.12) 

1, we can establish 



Q'"' : 



V\_ I'll + 



Pi 



(A.13) 



In the third line, we have defined the object Q^'^, which is neither a projector nor orthogonal 
to the -^ii^^'s but orthogonal to Pq'^. 

We are finally in a position to transform the polarization tensor for the parallel field 
configuration into its generalized form for arbitrary constant electromagnetic fields: 



w{k\A) = Uo pr + nil Pf + Pi"" + Q'^^ 

where no,||,_L and G are functions of the invariants and read: 



(A.14) 




2tt 



CO 










/ T 





-1 



-is0o 



eas ebs 



k'No 

sin eas sinh ebs ] N2'i^± — ^o'^y 



(A.15) 



Substituting the invariants into Eqs. (|A.4| ) and (|A.5|) , the functions Ni and 0o yield: 



No 

m 
N3 



m 



"^11 cosh ebs — cosh uebs ^ f ^ cos ueas — cos eas 



2 ebs sinh ebs 2 eas sin eas 

cosh uebs cos ueas — sinh uebs sin ueas cot eas cothe6s, 

cosh ebs — cosh uebs 
2 cos eas — 



2 cosh ebs 



sinh ebs 
cos ueas — cos eas 



sm eas 

1 — cos eas cos ueas cosh uebs cosh ebs — 1 



sm eas 



sinh efes 



+ sin ueas sinh i/e6s. (A. 16) 
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The scalars fjj are given by certain combinations of the invariants and can be found in 
Eq. (|A.10|) . The contact term given in Eq. ( |A.6| ) contributes equally to the Ilj's, 



c.t. = — e 



(A.17) 



but does not modify the function B, which is already finite. 

Note that Eq. ([A.14| ) almost appears in a diagonalized form except for the term B Q^^ . 



While Pq^ indeed projects onto an eigenspace of 11^'^ with eigenvalue Ho, this is generally 
not the case for the projectors -^['j^' QQ^'^. Although a further diagonalization is 

straightforward, we will not bother to write it down, since we only need the trace of 11^^ , 
which is simply given by: 



Ho + n,! + n 



Q% = 0. 



(A.18) 



In the actual two-loop calculation, the contact terms can be omitted for two reasons: 
first, it does not contribute to the thermal part, since the latter is finite; secondly, for 
the zero-temperature Lagrangian, a renormalization procedure is required anyway and, in 
particular, the mass renormalization would not be covered by an inclusion of the contact 
terms. 

Inserting Eq. ([A.15|) into Eq. (|A.18|) brings us to the explicit representation of the 
trace: 



a r ds r du e eas ebs 

^ 2ti J s J 2 a'^ + b'^ sin eas sinh ebs 

-1 



Zk{N2 - Ni) + e{2No{a^+b^) + b^N^ + a^N^) 

(A.19) 



This is the desired expression which is required in Eq. (0). For reasons of convenience, it 
is useful to rewrite the function (^q in terms of the variables k"^ and Zk- For this, we insert 
Eq. ( |A.10| ) into the first line of Eq. ( [A.16|) ; a reorganization yields: 



(A.20) 



where we implicitly defined: 

1 



A, 



IS 



2 a? + b'^ 
is 1 



2 a? + b'' 



cos ueas — cos eas cosh uebs — cosh ebs 
■ + ■ 



eas sm eas 



, cos veas — cos eas 



eas sm eas 



ebs sinh ebs 



, cosh uebs — cosh ebs 
ebs sinh ebs 



(A.21) 
(A.22) 



This provides us with the required necessities for the two-loop calculation in Sec. 2. 
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B Finite- Temperature Coordinate Frame 



In order to make the paper self-contained, we briefly review the construction of the finite- 
temperature coordinate frame as introduced in , and then apply it to the present problem. 

First, we define the vierbein e^^ which mediates between the given system labelled by 
yU, z/, ■ • ■ = 0,1,2,3 and the desired system labelled by the (Lorentz) indices A,B,--- = 
0,1, 2, 3 by: 



eo^ 



ei^ 



62^ 



where the quantity d abbreviates the combination of invariants: 

d := 2TS -g^ + 8^. 

The vierbein satisfies the identity 

(^b^ = 9ab = diag(-l, 1, 1, 1), 



(B.23) 



(B.24) 



(B.25) 



where qab ~ g^^ denotes the metric which raises and lowers capital indices. By a direct 
computation, we can transform the field strength tensors and the heat-bath vector: 



n 



Fab '■— ^AtiF^^eBu 



(T, 0,0,0), 

/ \f£ \ 

-v^ 

-^d/s -g/^/s 

V g/Vs / 



(B.26) 



Obviously, the new system corresponds to the heat-bath rest frame with the spatial axes 
oriented along the electromagnetic field in some sense. The components of the field strength 
tensor are now given by combinations of the invariants. 

In order to determine the form of Zk = k^F^^^k^F^ a = —k^FAcF'" sk^ , we need the 
square of the field strength tensor: 



'AB 



FacF'^b 



/ -S 

s- 
Vd 

V 



e 



VdG 







\ 



£ 



e 







£ 



-91 

£ 



(B.27) 
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This allows us to write Zk in the form: 



8 {k^f - 2Vdk'k' + {2T + E) ^[^-E] [k^ + 2^ k^k^ + ^ (P)^, 



(B.28) 



where fc^, /c^, k^ represent the components of the rotated momentum vector k^ = e^f^k'^. 

Now we can finally determine the desired form for the exponent in Eq. (|16]) in terms 
of finite-temperature coordinates: 



e 

(B.29) 



where again, /c", k^, k"^, k^ represent the components of k^ 



C Two-loop Effective Action of QED at Zero Tem- 
perature 

Dedicating this appendix to the derivation of the zero-temperature two-loop Lagrangian 
has two reasons: first, we want to make contact to well-known results, which serve as a 
check for our computations; secondly, our results will represent a generalization of the work 
of Dittrich and Renter 0] , who considered purely magnetic fields, to the case of constant 
electromagnetic fields. 

Here, we will only give a version of the unrenormalized effective Lagrangian, since the 
renormalization program requires detailed investigations which are beyond the scope of the 
present work. In particular, the mass renormalization has to be treated with great care 
|35| , taking a correct matching of the imaginary parts of the effective action into account. 

In Eq. ([2l|), we achieved a separation of the thermal and the zero-temperature parts 
in the integral Ii, concentrating on the zero-T case, we found: 

Il=' = TT^e-'-^^ (C.30) 



where qa = {Ak + a?Az) and qh = {Ak — b'^Az), and Ak and A^ are defined in Eq. (|A.21| ). 
We also need the second integral /|"^°, which is related to If^^ by Eq. ([T3|), leading us to: 



n2 c 



00 



'^^ - W J 





im's r „2 



a 



is' + qa^is' + qb) is' + qa){s' + qb)^ 



(C.31) 
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The integration over s' can be carried out with elementary techniques: 



rT=0 



167r2 



1 



Qaiqb - Qa) {qb - qaY Qa 



1 +. 1 ^^n*' 



Qbiqa - Qb) {qb - QaY Qa 

(C.32) 



Inserting these T = 0-contributions ( |C.3(]| ) and ( |C.32| ) into Eq. (|I^ and reorganizing the 
result a bit, we finally arrive at: 



oo 1 

a I ds r du 



(An) 



s J 2 

-1 



eas ebs 



sin eas sinh ebs 



X 



H : r : rrr m 



Qaiqb - Qa) QbiQa " Qb) (Qb " QaY Qa 



(C.33) 



where Nq, Ni, and N2 are functions of the integration variables and the invariants a and 
b, and are defined in Eq. (|A.16|) , whereas Qa and qi, after insertion of Eq. (|A.21|) can be 
written as: 



Qa 



IS cos ueas — cos eas 



eas sm eas 



Qb 



is cosh ebs — cosh uebs 
2 ebs sinh ebs 



{CM) 



Equation ( |C.33| ) represents our final result for the unrenormalized two-loop effective La- 
grangian of QED for arbitrary constant electromagnetic fields. In the limit of vanishing 
electric fields, we recover exactly the findings of 0; hence our comparably compact repre- 
sentation generalizes the results of |@] to arbitrary constant electromagnetic fields. 
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